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We investigate a class of one-dimensional, exactly solvable anisotropic XY spin-1/2 mod- 
' ' els in an alternating transverse magnetic field from an entanglement perspective. We find 

^ ' that a physically motivated Lie-algebraic generalized entanglement measure faithfully 

Q ^ portraits the static phase diagram - including second- and fourth-order quantum phase 

O ■ transitions belonging to distinct universality classes. In the simplest time-dependent sce- 

nario of a slow quench across a quantum critical point, we identify parameter regimes 
where entanglement exhibits universal dynamical scaling relative to the static limit. 

Keywords: Entanglement; Quantum Phase Transitions; Quantum Information Science 

' 1. Introduction 

Developing methodologies for probing, understanding, and controlling quantum 
' phases of matter under a broad range of equilibrium and non-equilibrium conditions 

. is a central goal of condensed-matter physics and quantum statistical mechanics. 

Since novel forms of matter tend to emerge in the deep quantum regime where ther- 
mal effects are frozen out, a key prerequisite is to obtain an accurate theoretical 
^ . , understanding of zero-temperature quantum phase transitions (QPTs).^ Aside from 

d ' its broad conceptual significance, such a need is heightened by the growing body 

of experimental work which is being performed at the interface between material 
science, quantum device technology, and experimental implementations of quantum 
information processing (QIP). Following the experimental realization of the Bose- 
Hubbard model in a confined ^^Rb Bose-Einstein condensate and the spectacular 
observation of the superfluid-to-Mott-insulator QPT,^ ultracold atoms are enabhng 
investigations into strongly interacting many-body systems with an unprecedented 
degree of control and flexibility - culminating in the observation of topological 
defects in a rapidly quenched spinor Bose-Einstein condensate.^ Remarkably, the 
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occurrence of a QPT influences physical properties well into the finite-temperature 
regime where real-world systems live, as vividly demonstrated by the measured 
low-temperature resistivity behavior in heavy- fermion compounds.^ 

From a theoretical standpoint, achieving as a complete and rigorous quantum- 
mechanical formulation as desired is hindered by the complexity of quantum corre- 
lations in many-body states and dynamical evolutions. Motivated by the fact that 
QIP science provides, first and foremost, an organizing framework for addressing 
and quantifying different aspects of "complexity" in quantum systems, it is natural 
to ask: Can QIP concepts and tools contribute to advance our understanding of 
many-body quantum systems? In recent years, entanglement theory has emerged as 
a powerful bridging testbed for tackling this broad question from an information- 
physics perspective. On one hand, entanglement is intimately tied to the inherent 
complexity of QIP, by constituting, in particular, a necessary resource for computa- 
tional speed-up in pure-state quantum algorithms.^ On the other hand, critically re- 
assessing traditional many-body settings in the light of entanglement theory has al- 
ready resulted in a number of conceptual, computational, and information-theoretic 
developments. Notable advances include efficient representations of quantum states 
based on so-called projected entangled pair states,^ improved rcnormalization-group 
methods for both static 2D and time-dependent ID lattice systems,^ as well as rig- 
orous results on the computational complexity of such methods and the solvability 
properties of a class of generalized mean-field Hamiltonians.* 

In this work, we focus on the problem of characterizing quantum critical models 
from a Generalized Entanglement (GE) perspective, by continuing our earlier 
exploration with a twofold objective in mind: first, to further test the usefulness 
of GE-based criticality indicators in characterizing static quantum phase diagrams 
with a higher degree of complexity than considered so far (in particular, multi- 
ple competing phases); second, to start analyzing time-dependent, non-equilibrium 
QPTs, for which a number of outstanding physics questions remain. In this context, 
special emphasis will be devoted to establish the emergence and validity of universal 
scaling laws for non- equilibrium observables. 

2. Generalized Entanglement in a Nutshell 
2.1. The need for GE 

Because a QPT is driven by a purely quantum change in the many-body ground- 
state correlations, the notion of entanglement appears naturally suited to probe 
quantum criticality from an information-theoretic standpoint: What is the struc- 
ture and role of entanglement near and across criticality? Can appropriate entan- 
glement measures detect and classify quantum critical points (QCPs) according to 
their universality properties? Extensive investigations have resulted in a number of 
suggestive results, see e.g. Ref. 11 for a recent review. In particular, pairwise entan- 
glement, quantified by so-called concurrence, has been found to develop distinctive 
singular behavior at criticality in the thermodynamic limit, universal scaling laws 
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being obeyed in both ID and 2D systems. Additionally, it has been established that 
the crossing of a QCP point is typically signaled by a logarithmic divergence of the 
entanglement entropy of a block of nearby particles, in agreement with predictions 
from conformal field theory. While this growing body of results well illustrates the 
usefulness of an entanglement-based view of quantum criticality, a general theoret- 
ical understanding is far from being reached. With a few exceptions, the existing 
entanglement studies have focused on analyzing how (i) bipartite quantum corre- 
lations (among two particles or two contiguous blocks) behave near and across a 
QCP under the assumption that the underlying microscopic degrees of freedom 
correspond to (ii) distinguishable subsystems (iii) at equilibrium. 

GE provides an entanglement framework which is uniquely positioned to over- 
come the above limitations, while still ensuring consistency with the standard 
"subsystem-based" entanglement theory in well-characterized limits. ^'^'''■^^ Physi- 
cally, GE rests on the idea that entanglement is an observer- dependent concept, 
whose properties are determined by the expectations values of a distinguished sub- 
space of observables f2, without reference to a preferred decomposition of the overall 
system into subsystems. The starting point is to generalize the observation that stan- 
dard entangled pure states of a composite quantum system look mixed relative to 
an "observer" whose knowledge is restricted to local expectation values. Consider, in 
the simplest case, two distinguishable spin-1/2 subsystems in a singlet (Bell) state, 

|ggii^ ^ I T)a'^U)b-U)a0| T)b ^ 

defined on a tensor-product state space H = Ha Hb- First, the statement 
that I Bell) is entangled - |Bell) cannot be expressed as \ijj)A <S^ \'p)b for arbitrary 
IV')^ € Ti-A, G ~ is unambiguously defined only after a preferred tensor- 

product decomposition of H is fixed: Should the latter change, so would entangle- 
ment in general. Second, the statement that |Bell) is entangled is equivalent to 
the property that (either) reduced subsystem state - as given by the partial trace op- 
eration, pA = TrB{|Bcll)(Bell|} - is mixed, Tr{p^} = 1/2(1 + Ea=..,j/,.(^a )') < 1' 
in terms of expectations of the Pauli spin-1/2 matrices acting on A. 

To the purposes of defining GE, the key step is to realize that a meaningful notion 
of a reduced state may be constructed for any pure state without invoking 

a partial trace, by specifying such a reduced "fi-state" as a list of expectations of 
operators in the preferred set O. The fact that the space of all 0-states is convex 
then motivates the following:^ 

Definition (Pure-state GE). A pure state \tp) e H is generalized unentangled 
relative to O if its reduced fl-state is pure, generalized entangled otherwise. 

For applications to quantum many-body theories, two major advantages emerge 
with respect to the standard entanglement definition: first, GE is directly appli- 
cable to both distinguishable and indistinguishable degrees of freedom, allowing 
to naturally incorporate quantum-statistical constraints; second, the property of 
a many-body state \ip) to be entangled or not is independent on both the choice 
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of "modes" (e.g. position, momentum, etc) and the operator language used to de- 
scribe the system (spins, fermions, bosons, etc) ~ depending only on the observables 
f2 which play a distinguished physical and/or operational role. 

2.2. GE by example 

For a large class of physical systems, the set of distinguished observables il may 
be identified with a Lie algebra consisting of Hermitian operators, 17 ~ f), which 
generate a corresponding distinguished unitary Lie group via exponentiation, I) 
G = e''' . While the assumption of a Lie-algebraic structure is not necessary for the 
GE framework to be applicable, ^'^^ it has the advantage of both suggesting simple 
GE measures and allowing a complete characterization of generalized unentangled 
states. In particular, a geometric measure of GE is given by the square length 
(according to the trace norm) of the projection of onto f): 

Definition (Relative purity). Let {Oe}, i — 1,...,M, be a Hermitian, or- 
thogonal basis for f), dim((]) = M. The purity of \ip) relative to f) is given by 



where K is a global normalization factor chosen so that < Pp, < 1. 

Notice that Pf, is, by construction, invariant under group transformations, that 
is, PijdV")) — P(j(G|?/')), for all G & G, as desirable on physical grounds. If, addi- 
tionally, f) is a semi-simple Lie algebra irreducibly represented on H, generalized 
unentangled states coincide^ with generalized coherent states (GCSs) of G, that is, 
they may be seen as "generalized displacements" of an appropriate reference state, 
|GCS({77£})) — exp{iJ2eVtOe)\ref). Physically, GCSs correspond to unique ground 
states of Hamiltonians in t): States of matter such as BCS superconductors or nor- 
mal Fermi liquids are typically described by GCSs. While we refer the reader to 
previous work^'^°'^'^ for additional background, we illustrate here the GE notion by 
example, focusing on two limiting situations of relevance to the present discussion. 

2.2.1. Example 1: Standard entanglement revisited 

The standard entanglement definition builds on the assumption of distinguishable 
quantum degrees of freedom, the prototypical QIP setting corresponding to N lo- 
cal parties separated in real space, and Ti = Hi (8) ... (8) Hn- Available means for 
manipulating and observing the system are then naturally restricted to arbitrary 
local transformations, which translates into identifying the Lie algebra of arbitrary 
local (traceless) observables, t)ioc = su(dim('Hi)) © . . . ©su(dim(HAr)), as the distin- 
guished algebra in the GE approach. If, for example, each of the factors TCi supports 
a spin-1/2, fijoc = spanja^ ]a = x,y, z,£ — 1, . . . , N}, and Eq. ([2]) yields 



M 




(2) 




(3) 
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which is nothing but the average (normahzed) subsystem purity. Thus, Ptji^^ quan- 
tifies multipartite subsystem entanglement in terms of the average bipartite entan- 
glement between each spin and the rest. Maximum local purity, P[, = 1, is attained 
if and only if the underlying state is a pure product state, that is, a GCS of the 
local unitary group Gioc — SU{2)i (g) . . . (g) SU{2)]y. 



2.2.2. Example 2: Fermionic GE 

Consider a system of indistinguishable spinless fermions able to occupy N modes, 
which could for instance correspond to distinct lattice sites or momentum modes, 
and are described by canonical fermionic operators Cj , on the 2^-dimensional Fock 
space Tipock- Although the standard definition of entanglement can be adapted to 
the distinguishable-subsystem structure associated with a given choice of modes (re- 
sulting in so-called "mode entanglement"), privileging a specific mode description 
need not be physically justified, especially in the presence of many-body interac- 
tions.^^ These difficulties are avoided in the GE approach by associating "gener- 
alized local" resources with number-preserving bilinear fermionic operators, which 
identifies the unitary Lie algebra u{N) — span{c]cj ; 1 < i, j < N} as the distin- 
guished observable algebra for fermionic GE. Upon re-expressing u{N) in terms of 
an orthogonal Hermitian basis of generators, Eq. ([2]) yields 



N , N 

4 



i<fe=l 3=1 

One may show^'' that a many-fermion pure state is generalized unentangled relative 
to u{N) if and only if it is a single Slater determinant (with any number of fermions), 
whereas -Pu(Af) < 1 fo'' ^^V/ state containing fermionic GE. Note that a Bell pure 
state as in Eq. ([T]) rewrites, via a Jordan- Wigner isomorphic mapping, in the form 

^ |T)a'8U)b-U)a'^|T)b ^ cl|vac)-4|vac) 
V2 V2 ' 

in terms of the fermionic vacuum |vac) = \1)a ^ U)s = U,i)- Thus, while |Bell) 
is maximally mode-entangled relative to the local spin algebra su(2) ® su(2), it is 
u(7V)-M7Tentangled - consistent with the fact that it is a one-particle state. 



3. Generalized Entanglement and Quantum Critical Phenomena 
3.1. Static QPTs 

Let us focus in what follows on a class of exactly solvable spin-1/2 one-dimensional 
models described by the following Hamiltonian: 



H 



— 2 — + — 2 — y y +2^ (-1) -^jc^z: (6) 



i=l 
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where periodic boundary conditions are assumed, that is, cr^ = cr^^^ ■ Here, 
7 € [0,1], h g [— cxDjCxd], and 5 € [— cxd,cxd] are the anisotropy in the XY plane, 
the uniform magnetic field strength, and the alternating magnetic field strength, 
respectively. For 6 = 0, the above Hamiltonian recovers the anisotropic XY model 
in a transverse field studied in Ref. 10, whereas (5 > 0, 7 = 1 corresponds to the 
Ising model in a alternating transverse field recently analyzed in Ref. 15. 

While full detail will be presented elsewhere, an exact solution for the energy 
spectrum of the above Hamiltonian may be obtained by generalizing the basic steps 
used in the standard Ising case,^'' in order to account for the existence of a two- 
site primitive cell introduced by the alternation. By first separately applying the 
Jordan- Wigner mapping to even and odd lattice sites, and then using a Fourier 
transformation to momentum space, Hamiltonian ((6]) may be rewritten as: 

H.^H..j:AlHa.. AV^{i| (1-^)}, 

k£K+ kGK+ 

where Hk is a four-dimensional Hermitian matrix, and = (a^, a_fe, is a 

vector operator, aj, denoting canonical fermionic operators that create a spinless 
fermion with momentum k for even (odd) sites, respectively. Thus, the problem 
reduces to diagonalizing each of matrices Hk, for k G K+. If efe^i, ek,2, efe,3, efe,4, with 
Cfc,! < efc,2 < < Cfc^a < tk.A denote the energy eigenvalues of Hk, then 

Hk = '^Knlln"fk,n , 

ri=l,...,4 

where 7^ jk,n are quasi-particle excitation operators for mode k in the nth band. 
At T = 0, the ek.i and ek,2 bands are occupied, whereas £^.3 and ek,4 are empty, 
thus the ground-state energy Egs = J2k£K+i''kA + £^,2), with ek,i < 0,£k,2 < 0. 

By denoting with |vac) the fermionic vacuum, and by exploiting the symmetry 
properties of the Hamiltonian, the many-body ground state may be expressed in 
the form |*)gs = UkeK+ 1*0, with 

(7) 

for complex coefficients determined by diagonahzing Hk, with J2l=i — ^■ 

Since QPTs are caused by non-analytical behavior of Egs, QCPs correspond to 
zeros of ek,2- The quantum phase boundaries are determined by the following pair 
of equations: = 6^ + 1; = +7^- The resulting anisotropic (7 > 0) quantum 
phase diagram is showed in Fig. [T] where, without loss of generality, we set 7 = 
0.5. Quantum phases corresponding to disordered (paramagnetic, PM) behavior, 
dimer order (DM), and ferromagnetic long-range order (FM) emerge as depicted. 
In the general case, the boundaries between FM and PM phases, as well as between 
FM and DM phases, are characterized by second-order broken-symmetry QPTs. 
Interestingly, however, Eqs develops weak singularities at 

(/i„ (5,) = (0, (5 = ±7) , iK,6c) = (±1,5 = 0), (8) 
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Fig. 1. Phase diagram of the spin-1/2 XY alternating Hamiltonian given in Eq. ((H) with 7 = 0.5. 

where fourth-order broken-symmetry QPTs occur along the paths approaching the 
QCPs (Fig.[Tl dashed-dotted hnes). In the isotropic hmit (7 = 0), an insulator- metal 
Lifshitz QPT occurs from a gapped to a gapless phase, with no broken-symmetry 
order parameter. For simplicity, we shall primarily focus on gapped quantum phases 
in what follows, thus 7 > 0. Standard finite-size scaling analysis reveals that new 
quantum critical behavior emerges in connection with the alternating fourth-order 
QCPs in Eq. Thus, in addition to the usual Ising universality class, charac- 

terized by critical exponents ly = l,z = 1, an alternating universality class occurs, 
with critical exponents — 2, z = 1. 

The key step toward applying GE as a QPT indicator is to identify a (Lie) 
algebra of observables whose expectations reflect the changes in the GS as a function 
of the control parameters. It is immediate to realize that Hamiltonian Eq. ([6]), 
once written in fermionic language, is an element of the Lie algebra so(2A'^), which 
includes arbitrary bilinear fermionic operators. As a result, the GS is always a GCS 
of so (2 A^), and GE relative to so(2A^) carries no information about QCPs. However, 
the GS becomes a GCS of the number-conserving sub-algebra u(iV) in both the fully 
PM and DM limit. This motivates the choice of the fermionic u(-/V)-algcbra discussed 
in Example ^ as a natural candidate for this class of systems. Taking advantage of 
the symmetries of this Hamiltonian, the fermionic purity given in Eq. ([4]) becomes: 



Analytical results for Pu(Ar) are only available for (5 = 0, where GE sharply detects 
the PM-FM QPT in the XY model. 1° Remarkably, ground-state fermionic GE still 
faithfully portraits the full quantum phase diagram with alternation. First, deriva- 
tives of Pu(7v) develop singular behavior only at QCPs, see Fig. [2] (left). Furthermore, 




keK+ 



(9) 



} 
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GE exhibits the correct scaHng properties near QCPs.^° By taking a Taylor expan- 
sion, Pu{N) (h) — Pu(N){hc) ~ ~ {h — hcY, where ^ is the correlation length, the 
static critical exponent v may be extracted from a log-log plot of Pu{n) for both 
the Ising and the alternating universality class, as demonstrated in Fig. [5] (right). 




static field strength, h In |h-h^| 



Fig. 2. Pii(N) 3s ^ static QPT indicator. Left panel: Purity and rescaled purity derivative vs 
magnetic field strength. Inset: second derivative for N = 1000,2000,4000,8000 (top to bottom). 
Right panel: Determination of i/ for both the alternating and Ising (inset) universality class. 



3.2. Dynamic QPTs 

While the above studies provide a satisfactory understanding of static quantum crit- 
ical properties, dynamical aspects of QPTs present a wealth of additional challenges. 
To what extent can non-equilibrium properties be predicted by using equilibrium 
critical exponents? The simplest dynamical scenario one may envision arises when a 
single control parameter is slowly changed in time with constant speed Tq > 0, that 
is, g{t) — gc = [t — tc)lTq, so that a QCP is crossed at t = tc (ic = without loss 
of generality). The typical time scale characterizing the response of the system is 
the relaxation time t = h/ A ^ \g{t) ~ 9c\~'"^, A being the gap between the ground 
state and first accessible excited state and z the dynamic critical exponent.^ Since 
the gap closes at QCPs in the thermodynamic limit, r diverges even for an arbi- 
trarily slow quench, resulting in a critical slowing- down. According to the so-called 
Kibble-Zurek mechanism (KZM),^^ a crossover between an (approximately) adia- 
batic regime to an (approximately) impulse regime occurs at a freeze-out time — 
whereby the system's instantaneous relaxation time matches the transition rate, 

T{i) = \{g{i)-g,)/g'{i)\, t-r-A^+D, 

resulting in a predicted scaling of the final density of excitations as 



(10) 
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While agreement with the above prediction has been verified for different quantum 
systems/^ several key points remain to be addressed: What are the required phys- 
ical ingredients for the KZM to hold? What features of the initial (final) quantum 
phase are relevant? How does dynamical scaling reflect into entanglement and other 
observable properties? 

In our model, the time-evolved many-body state at instant time t, = 
nfeei<-+ \^k{t)), may still be expressed in the form of Eq. ([7]) for time-dependent 
coefficients a — 1,...,6, computed from the solution of the Schrodinger 

equation, subject to the initial condition that \^(t — oo)) = |\E'g5(— oo)). The fi- 
nal excitation density is then obtained from the expectation value of the appropriate 
quasi-particlc number operator over the final state, 

n{tF) = ^mF)\ J2 (7l,37fc,3 + 7i,47fc,4) mtp)) ■ 

keK+ 

As shown in Fig. [3] (left), the resulting value agrees with Eq. dTUI) over an appropriate 
Tg-range irrespective of the details of the QCP and the initial (final) quantum phase: 

nitpY''"^ - , „(i^)^ite™at*«g _ -2/3 




Fig. 3. Dynamical scaling of the excitation density. Left panel: log-log plot for Ising universality 
class (FM to PM). Right panel: alternating universality class (FM to FM), with log-log scaling 
plot in the inset. 

Remarkably, however, our results indicate that scaling behavior holds throughout 
the entire time evolution (see Fig. [3l right), implying the possibility to express the 
time-dependent excitation density as: 

„(t)=r-/(-+i)F(^), 

where is a universal scaling function. Numerical results support the conjecture 
that similar universal dynamical scaling holds for arbitrary observables.^^ In par- 
ticular, fermionic GE obeys scaling behavior across the entire dynamics provided 
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Fig. 4. Dynamical scaling of Pi,(A^) for the alternating and the Ising (inset) universality class. 

that the amount relative to the instantaneous ground state \^{t))Qs is considered: 
AP,(^.)(t) ^ Puw(|$W)) - Puw(|*(t)>Gs) - r-''/(-+i)G(^) , 

for an appropriate scaling function G, see Fig. 21 

It is important to stress that the above discussion apphes to control paths which 
originate and end in gapped phases. In the isotropic limit 7 = 0, we observe no 
scaling of the form Eq. pO|) if the system is driven to/from the superfluid gapless 
phase. 

4. Conclusion 

In addition to further demonstrating the usefulness of the GE notion toward char- 
acterizing static quantum critical phenomena, we have tackled the study of time- 
dependent QPTs in a simple yet illustrative scenario. Our analysis points to the 
emergence of suggestive physical behavior and a number of questions which deserve 
to be further explored. In particular, while for gapped systems as considered here, 
the origin of the observed universal dynamical scaling is likely to be rooted in the 
existence of a well-defined adiabatic (though non-analytic) limit - as independently 
investigated in Ref. 20, a rigorous understanding remains to be developed. We ex- 
pect that a GE-based perspective will continue to prove valuable to gain additional 
insight in quantum-critical physics. 
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